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Summary 

Theoretically, the electric field becomes infinite at corners of two and three di- 
mensions and edges of three dimensions. Conventional finite-element and bound- 
ary element methods do not yield satisfactory results at close proximity to these 
_ ^ singular locations. In this paper, we describe the application of a fast and accurate 

Q ■ BEM solver (which uses exact analytic expressions to compute the effect of source 

tyj , distributions on flat surfaces) to compute the electric field near three-dimensional 

. y I corners and edges. Results have been obtained for distances as close as Ifim near 

^' the comer / edge and good agreement has been observed between the present results 

(— I . and existing analytical solutions. 
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Introduction 



►^ I Accurate computation of electric field near comers and edges is important in 

Cf^ • many applications related to science and engineering. While the generic problem 

^ . is important even in branches like fluid and solid mechanics, the specific electro- 

^ I static problem is very important in the modeling and simulation of micro electro- 

O ' mechanical systems (MEMS), electric probes and antennas, to name a few. While 

it is true that the ideal corner / edge does not exist in nature, the singularity being 
smoothed by rounding corners, sharp increase in charge density does occur near 
Q ' these geometries. Since the electric field is proportional to the square of the charge 

* ^ ■ density, it is very important to estimate the charge density and the resulting varia- 

^! tion of electric field in their vicinity. 

^ , 

"' There have been many attempts at solving the above problem using the finite- 

element (FEM) and the boundary element method (BEM). However, because of 
^ . the nature of the problem, significant modifications to the basic method needed to 

H I be carried out to handle the boundary singularities. On the simpler side, special 

mesh refinement schemes have been used while on the more sophisticated side, 
the form of local asymptotic expansions (which may often be found) have been 
used. Very effective FEM solvers have been developed Q which calculate singular 
coefficients by post-processing the numerical solution. These solvers improve the 
solution by refining the mesh and changing the degree of piecewise polynomials. 
Similarly, several accurate BEM solvers have been developed Q, ||3l- For exam- 
ple, in the singular function boundary integral method, the singular coefficients are 
calculated directly and the solution is approximated by the leading terms of the 
local asymptotic solution. 
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In this paper, we present a solution to the comer / edge problem by using a re- 
cently developed three-dimensional BEM solver ||4]|. This solver uses exact analytic 
expressions for computing the influence of charge evenly distributed on discretized 
flat elements on the boundary. Through the use of these closed-form expressions, 
the solver avoids one of the most important approximations of the BEM, namely, 
the assumption that the charge on a boundary element is concentrated at one point 
(normally the centroid) of the element. As a result, the computed capacitance ma- 
trix is very accurate and the solver is able to handle difficult geometries and charge 
distributions with relative ease. The solver has been used to solve the classic prob- 
lem of two planes intersecting at various angles. Exact solution to this problem 
exists 151 and our results have been compared with the exact results. The compar- 
ison is very good even as close as 1 /um to the comer / edge. Especially important 
is the fact that the solver produces quite accurate results even for the case of an 
edge. It has also been possible to reproduce the dependence of the strength of the 
electric field as a function of the distance from the geometric singularity. All the 
calculations have been carried out in three dimensions and, here, we also present 
the variation of the electric field along the corner or edge. It is observed that the 
electric field retains its mid-plane value for much of the distance along its length, 
but increases significantly within 20% of the axial ends. It may be noted here that 
for these calculations, only algebraic mesh refinement near the edge was applied 
and no other sophisticated techniques such as those mentioned above were applied. 
This fact made the development of the solver and its application free of mathemat- 
ical and computational complications. Since corners and edges play an important 
role in many electro-mechanical structures, the solver can be used to study the elec- 
trostatic properties of such geometries. 

Theory 

We have considered the geometry as presented in (51 in which two conducting 
planes intersect each other at an angle p. The planes are assumed to be held at a 
given potential. A circular cylinder is also included that just encloses the two in- 
tersecting plane, has its center at the intersection point and is held at zero potential. 
The general solution in the polar coordinate system (p, (|)) for the potential (<I>) close 
to the origin in this problem has been shown to be 



^(pA) = ^ + L a^p""^/^m(m7l(^/P) (1) 

m=l 

where the coefficients a^ depend on the potential remote from the corner at p = 
and V represents the boundary condition for <I> for all p > when ([) = and (|) = 
p. In the present case where a circular cylinder just encloses the two plates, the 
problem of finding out a,„ reduces to a basic fourier series problem with a well 



known solution 



a,n = — for odd m (2) 

mn 

It may be noted here that the series in Q involves positive powers of p'^/P, and, 
thus, close to the origin (i.e., for small p), only the first term in the series will be 
important. The electric field components (Ep , E,^) are 



Ep{p,<^) = -? L a,p("™/W-i««(m7l(^/P) (3) 



E^{p,<\,) = -^ £ «,„p('"'^/P'-icos(m7l^/P) (4) 

r^ m= 1 



Thus, the field components near p = vary with distance as p('^/P) ^ and this fact 
is expected to be reflected in a correct numerical solution as well. 

Numerical Solution 

While the above theoretical solution is a two-dimensional one, we have used 
the BEM solver to compute a three-dimensional version of the above problem. In 
order to reproduce the two-dimensional behavior at the mid-plane, we have made 
the axial length of the system sufficiently long, viz., 10 times the radius of the 
cylinder. The radius of the cylinder has been fixed at one meter, while the length 
of the intersecting flat plates has been made a micron shorter than the radius. The 
length of the plate has been kept smaller than the radius of the cylinder to avoid the 
absurd situation of having two values of the voltage at the same spatial point. We 
believe that it has been possible to maintain such a small gap as 1 /um between the 
circular cylinder and the flat plates without facing numerical problems because the 
BEM solver we have used computes the capacity matrix using analytic expressions 
which calculate accurate values of potential and electric field at extremely close 
distances from the singular surfaces IH. 

The cylinder has been discretized uniformly in the angular and axial directions. 
The flat plates have also been uniformly discretized in the axial direction. In the ra- 
dial direction, however, the flat plate elements have been made successively smaller 
towards the edges using a simple algebraic equation in order to take care of the fact 
that the surface charge density increases considerably near the edges. The elec- 
tric field has been computed for various values of P (1.2571, 1.571, USnandlii), 
referred to as cases 1,2,3 and 4 respectively in the following section. 



~C3 
CD 



1000 



100 



10 



0.1 





□ 


■ 1 


■ 1 


' ■ 1 


case 1 


+ 


: 




^^ 








case 2 


^ 


: 






" 's,^. 






case 3 




- 


r- 


X . 


^< 






case 4 


□ 


-. 








"Q, 








: 






X -. 


^ "■- 








: 




+ ---_^ 




X -. 








- 


~ 




~~~~4---- 


~~~~H-~^^^ 


X-. 


■ --- "B 




~ 










~~^-^~ 


^^_ -M: -., 




- 












~~"t~--:"; 


^*,__ 


,j 










. . 1 




. . 1 


'- 



1e-07 1e-06 1 e-05 1 e-04 0.001 0.01 0.1 
Distance from corner (m) 

Figure 1 : Comparison of electric field values 



Results 

In Figur^ we have presented a comparison of the variation of the electric field 
as a function of the distance from the edge as found from the analytical solution ^ 
and the BEM solver. Computations have been carried out upto a distance of 1 /um 
from the edge and to properly represent the sharp increase in field, logarithmic 
scales have been used. The computed electric field is found to be in remarkable 
agreement for all values of p. There is a small disagreement between the two 
results only at the point closest to the corner / edge (« 1.3%). At present, we are in 
the process of improving the BEM solver so that this error can be minimized. 

Although, the distance dependence of the electric field is expected to match 
with the theoretical prediction because the computed and the theoretical estimates 
are found to match closely, in Figure |2 we have plotted curves corresponding to 
pJt/p-i aj^(j compared the computed values as points against the curves. Thus, 
we can confirm that the computed electric field obeys the p('^/P)^^ relation quite 
accurately near the corner / edge. Here, to emphasize and visualize the sharpness 
of the increase of electric field near the singular geometry, only the distance has 
been drawn using the logscale. 

In Figur^U we have presented the variation of the electric field at various dis- 
tances from the corner / edge along the length of the geometry. While the field is 
found to increase considerably towards the end, for most of the length (^ 80%), the 
two-dimensional value seems to be a good approximation for the real value. How- 
ever, the deviation (as large as ^ 50%) from the analytic solution towards the front 
and back ends of the geometry should be taken into consideration while designing 
real-life applications. At close distances near the singular location, oscillation in 
the electric field value is observed. Initial study indicates that the source of this 
oscillation is possibly the large size of the boundary elements in the axial direction. 
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Figure 2: Dependence of electric field as a function of distance 
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Figure 3: Variation of electric field along the axis 



Distance 


Exact 


Computed 


Error (%) 


0.1 


1.83015317 


1.829997 


0.0085332 


0.01 


6.303166063 


6.302387 


0.012359868 


0.001 


20.11157327 


20.103534 


0.039973352 


0.0001 


63.65561168 


63.548305 


0.168573794 


0.00001 


201.31488353 


199.998096 


0.654093481 


0.000001 


636.6191357 


628.098310 


1.33844951 



Table 1 : Comparison of electric field values 

In order to facilitate numerical comparison, we have presented the values of 
the electric field in the mid-plane as obtained from the exact solution and the BEM 
solver in Table[2for the three dimensional edge. We have chosen this case since it 
is known to be the most difficult among all those considered in this work. It is clear 
from the table that the present solver does indeed yield very accurate results. 

Conclusions 

A fast and accurate BEM solver has been used to solve the complex problem of 
finding the electrostatic force field for three dimensional corners and edges. Accu- 
rate solutions have been obtained upto distances very close singular geometry. The 
two dimensional analytic solution has been found to be valid for a large part of the 
geometry, but near the axial ends, the difference has been found to be significant. 
The results and the BEM solver is expected to be very useful in solving critical 
problems associated with the design of MEMS, integrated circuits, electric probes 
etc. Since these problems are generic in nature, the solver should be important for 
analysis of problems related to other fields of science and technology. 
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